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' Abstract 

X) ' 

, In this paper we study the controUabihty of an artificial advection-diffusion system through the 

' boundary. Suitable Carleman estimates give us the observability on the adjoint system in the one 

dimensional case. We also study some basic properties of our problem such as backward uniqueness 
and we get an intuitive result on the control cost for vanishing viscosity. 
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O ■ Introduction 

o 

Artificial advection-diffusion problem In the present paper we deal with some advection-diffusion 
problem with small viscosity truncated in one space direction. Our interest for the linear advection 
diffusion equation comes from the Navier-Stokes equation, but it arises also in other fields as, for example, 
meteorology. For a given viscosity e > 0, the incompressible Navier-Stokes equation can be written as 

ft + {f-V)f - eA/ + Vp - 
divif) = 

OO . where / is the velocity vector field, p the pressure, V the gradient and A the usual Laplacian. Considering 

' the flow around a body, we have that / is almost constant far away from the body and equal to a (see 

. |10|). Linearizing the equation, we get the following equation for the vorticity 

(N 

O I ut + a.Vu - eAu = 0. 

ON 

. . , In the sequel, we assume for simplicity that a is the nth unit vector of the canonical basis of ii^ . 

_ ^ ' When one computes the solution of this problem, one can only solve numerically this problem on a 

, bounded domain. A good way to approximate the solution on the whole space may be given by the use of 

' artificial boundary conditions (see For any T > 0, we hence consider the following control problem 

. 5^ , on = M"^^ X (— L,0) {L some positive constant) 



ut + dnii - eAu = in (0, T) x fi, 

e{ut + d^u) = V on (0, T) x Fq, 

e{ut + d^u) + u = on(0,r)xFi, 

u(0, .) = uq in 51, 



where Fq := M"^^ x {0} and Fi := R"^^ x {—L} forms a partition of the boundary dil. Here we have 
denoted 9„ the partial derivative with respect to a;„ and di, the normal derivative. 
We are interested in the so-called null controllability of this system 

for given mq, find u such that the solution of (5t,) satisfies u(T) = 0. 

Using classical duality arguments, we will be interested on the observability of the adjoint system through 
To. 
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Let us first introduce the adjoint system: 




(ft + dnip + eAip = in (0, T) x ft, 
e{^t ~ d,^) - ^ = on(0,r)xro, 

^t-a^v^^o on(o,r)xri. 



ip(T, .) = (fiT in ^■ 



The observabihty inequahty corresponding to the previous controUabihty property is: 



there exists C > such that ||(^(0, .)\\x < C'||(/3||i2((o.T)xro) ^Vt G X, 



(1) 



where the space X will be defined below. 

One can notice that, when the viscosity e vanishes, our heat system {Sy) tends to a transport system 
with Dirichlet boundary condition on Fi. This phenomenon of degeneration of a parabolic problem to an 
hyperbolic one has been studied in several papers: see, for instance, [2] (one dimensional heat equation) 
and [5] (Burgers equation) . Similar results of interests can also be found in [1] . 



We will denote by Cotsi^) the cost of the null-control, which is the smallest constant C which fulfills 
the observability estimate ([T]). Our main result of the paper is the following: 

Theorem 1 Assume n ~ 1 and T,L > 0. 

• There exists C > such that the solution of problem (S") satisfies (Qp. Consequently, for every 
Uq £ X , there exists a control v G L^{{0, T) x Tq) with 



such that the solution of the problem (Sy) satisfies u{T) = 0. 

• Furthermore, ifT/L — 1 is large enough, the cost of the null-control Cotsi^) tends to zero exponen- 
tially as s ^ 0. 

Remark 1 One can in fact obtain an observability result for the adjoint system on Ti, that is 

\\ip{0,.)\\x < C||<^|U2(r„x(o,T)), 

for any ip solution o/(S"). This provides some controllability result for the direct system on Ti: one can 
find a function v depending continuously on uq so that the solution of 



satisfies u{T, .) = 0. 

Remark 2 The fact that the control cost vanishes tells intuitively that the state is almost null for T/L — 1 
big enough. This is to be connected with the fact that, for e = 0, the system is purely advective and then 
that, for T > L, its state is null. 

In some context of inverse problems (be able to know the origin of a polluted river for instance), it 
can be interesting to know if the observation of the solution of the direct problem on the boundary part 
Fi or Fq can allow us to recover the initial data. The corresponding result is presented now and will be 
proved at the end of the first section. 

Proposition 2 If n — \ and T > 0, if the solution of (Sq) with initial data Uq (z X satisfies u = on 
(0, T) X Ti, then uq = 0. However, there is no constant C > such that the following estimate holds true 



Main results We define X as the closure of T>{^) for the norm \\u\\x ■ 




M|L2((0,T)xro) < C'll"o||jf 



ut + dxU — edxxU = in (0, T) x $7, 
ut + dyu = on (0, T) x Fq 

e{ut + dyii) + u = v on {0,T) X Fi 




|uo|U < C||M||L2((o^T)xri) Vmo e X. 



(2) 



3 



Remark 3 Using Remark 1, one can also obtain similar result for Tq. 

The rest of the article is organized as follows: in the first section, we show the well-posedness of 
the direct and the adjoint problems using some semi-groups approach. In the second section, we will 
adapt Carleman inequality to the case of our one-dimensional problem. The third section is intended to 
explain how to get observability in the onc-dimcnsional case and the equivalence between observability 
and control. 

Notation: 

A < B means that, for some universal constant c > Q, A < cB. 

A ^ B means that, for some universal constant c > 1, c~^B < A < cB. 



1 Well-posedness and basic properties of systems 
1.1 Homogeneous problems 

We will use some semi-group results to show existence and uniqueness of the homogeneous direct problem 
(that is ("S'o)). This will enable to define solutions of the system (Sq) as a semigroup value. We define 
H = {X, ||.||a'),^ = H^{^) endowed with the usual norm 11. 1 1 and we consider the bilinear form on V 
defined by 



a(Mi,U2)=e / Vui.Vit2 + / 0,iMiM2 + / itiW2- (3) 
Jn Jn Jti 

Using Riesz representation theorem, one can define an operator A such that a(wi,U2) =< — ^wi,U2 >x 
and 'D{A) = {u £ X; — eAu -I- 9„u £ L'^iP), d^u £ L?{dil)} equipped with the graph norm 

II II II II \o-{u,v)\ 
\\u\\v(A) ^ \\u\\x + sup —— . 

(So) might now be written in the following abstract way 

ii = Au in X, 
w(0, .) = uq in X. 

Proposition 3 A generates an analytic semi-group (e*'^)t>o on X . 

Proof Using Hillc-Yoshida theorem, we will show that ^ is a maximal monotone operator. 
• First, Green- Riemann formula gives 

< Au,u>x= -e [ \\7u\'^ -I- [ lup, 
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SO the monotony is proved. 

Given v G X, we have to solve (/ — A)u = v, that is to find u e T>{A) such that 

Vu' e X,< {I - A)u, u' >x=< V, u' >x ■ 

The left-hand side term of this equation is a continuous bilinear form B on the space V while 
the right-hand side term is a continuous linear form L on V. Moreover, using ([3]), one can easily 
compute 

BM = \\u\\l+e I |Vup + i / |«p>min{l,£}hf, 
Jn ^ Jdn 

which show that B is coercive on V. Consequently, Lax-Milgram theorem shows that there exists 
u E V such that B{u,u') = L{u'), Vu' £ V. Using test functions u' £ T>{Q) additionnaly gives 
u £ ^{A) and our proof ends. 
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In particular, for every initial data uq G X, we have existence and uniqueness of a solution u G C(R+, X) 
to (5*0)- We will call these solutions weak solutions opposed to strong solutions i.e. such that uq G 2?(-4) 
and which fulfill u G C(M+,X»(y4)) n C^(R+,X). One can notice that, using the density of ViA) in X, a 
weak solution can always be approximated by a strong solution. 

Concerning the dual system, one can show that it may be written in the following abstract way 

( (p + A*ip = in X, 
\ ip{T, .) = LfT in X. 

In the same way as above, the adjoint operator A* is shown to be maximal monotone with domain 'D{A*). 
Thus, for every initial data ipx G X, we have existence and uniqueness of a solution tp G C(M.'^,X) to 
(5") given by means of the backward semigroup (e^^"*-*"^ )t>o- We will also speak of weak solutions or 
strong solutions in this situation. 



1.2 Nonhomogeneous direct problems 

We consider a slightly more general system 



go.gi, 



Ut + dnU — eAu ~ f in (0, T) x VI, 

e{ut + d^u) = go on (0, T) x Tq, 

e{ut + d^u) + u~gi on (0, T) x Fi, 

w(0, .) — Uq in n, 



with / G L^{{0, T) X n), go G L^{{0, T) x Tq) and 51 G ^^((o, T) x Ti). We consider a function g on dVL 
such that g = gi on T i (i = 0, 1). 

Definition 4 We say that u G C([0,r],X) is a solution of (Sf^gg^g-^) if, for every function 
ip G C([0,r] ,ff2(r2)) nCi([0,T] ,X), the following identity holds true 

u{ipt + dn(p + sAip) + / / u{e{ipt ~ d„ip) - lr„y:>) + / f^+ I 9V 



Jn Jo Jdn Jo Jn Jo Jen 



u{T)ip{T) - uoifiO) + e u{T)ip{T)~e uoip{0). 
n Jn Jan Jan 

Proposition 5 Let T > 0, uq e X , f e L'^{{0,T) x n), go G L^((0,T) x Tq) and gi G i^((0,r) x Ti). 
Then (S'/.g^^gj) possesses a unique solution u. Moreover, one has the following estimate 

I"I|l°°((o,t),x) < C{T,e) (||fiolU + ll/IU2(ox(o,T)) + \\9\\L^{dnx{o,T))) , 
with C{T,e) only depending on T and e. 

Proof Using Riesz representation theorem, we first define F{t) & X such that 



<F{t),u>x= / f{t)u+ / g{t)u, WueX. 
Jn Jdn 

It is now easy to check that u is solution to (S'/.go,gi) if and only if 

u{t) = e*^uo + f e'^'"^-^F{s)ds, Vt G [0, T]. 



Finally, using classical estimates of the norm of e*"^, this expression is well-defined and we obtain the 
required estimate. □ 

We now focus on some regularization effect for our general system and for some technical reasons we 
want to have some explicit dependence of the bounds on e. The result is given below. 
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Lemma 6 Let e e (0,1), / G L^{{0,T) x n), go G L^i{0,T) x Tq) and gi G L^iiO,T) x Ti). 

• Let uq ^ X . Then if u is the solution of {S f,gg^g-^), u belongs to L^{{Q,T), H^{il)) nC{[Q,T], X) and 

£^^'^\\u\\Lmo,T)MHn)) + I|w||l-((o,t),x) < C'(I|mo||x + ||/||L2((o,T)xn) + llffllL2((o,T)xan)) (4) 
for some C > on/?/ depending on T . 

• Lei now Mo e T'/ie" « e L2((0, T), n C([0, T], e ^^(-(-o^ 7^)^ x) and 

£^/^||u||ioo([o,T].ffi(n)) + £(I|Am||i,2((o_t)xl2(o)) + ||9^u||L2((o^T)xan))) + ||9tu||i2([o.T],x) 

C (5) 

- ^T72'^ll"oll-H"'(f^) + ll/IU2((o,T)xa) + ll.9llL2((o,T)xaji)), 

where C > onZj/ depends on T . 
Proof 

• First, we multiply our main equation by u and we integrate on fi. An application of Green- Ricmann 
formula with use of the boundary conditions gives the identity 

\dt\\u\\\ + \ [ \u\^ + e I iVnp ^ [ fu + [ gu^ 



2 " 2 
which immediately yields 



an Jn Jn Jon 



dthVx + J^^^ W\' + 2e 1^ |Vu|^ < \\u\\j, + + Mi.^eny 



Finally, Gronwall's lemma yields ([4]). 



• Now, we multiply the main equation by dtu and we integrate on which, using Green-Riemann 
formula and the boundary conditions provide the identity 

[2 I / Q „,C) „, I ^ Q / / |V7„,|2\ I 1 Q / / L,|2 



\\dtu\\j,+ / dnudtu+-dt[ / \Vu\'\+-dt[ / = / fdtu+ / g^tu, 

'o ^ \Jo / ^ \Jo / Jn Jon 



Now 



•J il, «^ O 

and, using again Young's inequality, we have 

fdtu+ [ gdtu^\\\dtu\\l + Cml^n) + \\9\\h(en))- 
We have thus proven the following inequality 

119*^111 + 2sdt (^^^ |Vnp^ + 2dt ^ |np < 2 J^^ \Vu\' + C (||.9lli.(a^^) + WfWl^n)) ■ (6) 

Gronwall's lemma and (jH) provide 



^^1' =^ § ^ + ll.9Wlli2(ao))rfi+ ll"o||^i(^j)J , (7) 

for t G (0, r). We now inject this into ^ to get the second part of the required result 



T 



T 



\\dtu{t)\\\dt I iWimhin) + \\9mhidn))dt+\\uo\\'jii^n) 
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Writing the problem as 

- Au = i (/ - dtu - dnu) in (0, T) x fi, 

dyu=^-dtu on(0,T)xro, 

dyu^^-dtu~\u on(0,r)xri, 

we get that u e i^((0, T), 2?(^)) and the existence of some constant C > such that 

ll^"llL2((0,T)xf2) + ll5,yU||L2((o, T)xaa) 

c 

< — (II/I|l2((o,T)xO) + ||.9lU2((o,T)xaO) + ll'^l|L2((0,T),ffi(n)) + \\dtu\\L'i({0,T).X)) 

which gives the last part of the required result ([5]). 



Remark 4 Working on the non-homogeneous adjoint problem and using the backward semigroup 
^(T~t)A ^ pj^g ij^^g show that the estimates of Lemma\^hold true for the system 

{(pt + dn(p + sA(p = f in{0,T)xn, 
e{tpt ~ d^ip) -tp^ go on (0, T) x Tq, 
e(v?t - a,</j) = <7i on(0,T)xri, 
f{T, .) = </?T in 

1.3 Backward uniqueness 

We will show here the backward uniqueness of systems (5*0) and (S"), thanks to the well-known result of 
Lions-Malgrange (|Hj) and the regularization effect. 

Lemma 7 Let uq ^ X and S G (0,T). Then, the solution u of (Sq) satisfies u G C{[5,T]^'D{A) . 

Proof Our strategy is to show that u G L^{{t,T); H'^{n)) and dtu G L'^{{t,T); H'^{n)). This easily 
implies that u G C{[t,T]; H^{n)). 

Wc first select some regular cut-off function 61 such that 6*1 = 1 on {6/2, T) and 9i = on (0, S/i). If 
ui = 9iu, we have that 

dtUi + dnUi — eAui ~ 9[u in (0, T) x fl, 
s{dtui + di,ui) — s9[u on (0, T) x Tq, 

s{dtui + di,ui) + ui — e9[u on(0,T)xri, 
Mi(0, .) = inn, 

that is ui satisfies {Sg'^u,E9{u,ee{u)- An application of Lemma [S] gives that ui G L'^{{0,T),'D{A}) and 
dtUi G L^{{0,T),X). This implies in particular that u G ^^((J^ t), 

Wc now focus on dtU and we select another cut-off function 6*2 such that 6*2 = 1 on (5, T) and 6*2 = on 
(0,(5/2). We write the system satisfied by 92U and we differentiate it with respect to time. One deduces 
that U2 = dt{92u) is the solution of 

{Sgi^u' +e'^ii,e(e'^u' +0'^u),e(0'^u' +e'^u)) ■ 

Using that 9^ is estimated by 9i, LemmaHgives that U2 G ^^((0, T), V{A)) that is dtU G L'^{{5, T), V{A)). 

□ 



Proposition 8 Assume that u is a weak solution of (So) such that u(T) = 0. Then uq = 0. 
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Proof li < S < T, then Lemma [7] shows that u{S) G ^{A). We will now apply Theorcme 1.1 of [S] to 
u as a solution of (So) in the time interval {S, T) to show that .) = 0. The bilinear form a defined in 
can be split into two bilinear forms on V defined by 

a(ui,U2) = ao(t, ui,U2) + ai{t,ui,U2) 

with 



ao(i,Ml,U2) = £ / VU1.VU2, ai(t,Ml,U2)= / 9„UiM2 + / UlU2. 

The four first hypotheses in [S] (see (1.1)-(1.4) in that reference) are satisfied since u G L^{6,T;V) fl 
H^{S,T;H) (from Lemma |6] above) , u{t) e 2?(yl) for almost every t € ((5,T), wt = Au and u(r, .) = 0. 
On the other hand, it is clear that ag and ai are continuous bilinear forms on V and that they do not 
depend on time t. It is also straightforward to see that 

V-u G y, ao{t,u,u) + ||u||^ ^ min{l,£}||u|p 

and, for some constant C > 0, 

yu,v G V, \ai{t,ui,U2)\ C||ui||||u2||x. 

This means that Hypothesis I in that reference is fulfilled. We have shown that u{d, .) = 0, which finishes 
the proof since 



uq ~ lim u(5) = 0. 



□ 



Remark 5 The result also holds for the adjoint system: ifip is a weak solution o/(S") such that ip{0) = 
then ifT = 0. The proof is very similar to the one above and left to the reader. 

1.4 Proof of Proposition [2] 

The fact that u = on (0, T) x Fi implies uq = is a straightforward consequence of the backward 
uniqueness of system (Sq)- Indeed, similarly as ([T]), one can prove the observability inequality 

\HT,.)\\x < \\u\\L^{0,T)x^^), 

(see Remark [6] below), which combined with Proposition [8] gives uq = 0. 

To show that ^ is false, we first need to prove some well-posedness of (Sq) with uq in a less regular 
space. For s g]0, l/2[, we define X'^ as the closure of I?(f2) for the norm 



where H'^{^1) (resp. H''{dil)) stand for the usual Sobolev space on D, (resp. dft). One easily shows that 
X" is a Hilbert space. We now denote the set of functions U2 such that the linear form 

ui G T>{ri) I — >< ui,U2 >x 

can be extended in a continuous way on X'^. If ui G X~'^, U2 G X" we denote this extension by 

< Ml,U2 >-s, fi- 
ll uo G X^^, we say that u is a solution by transposition of (So) if, for every / G L^((0,T) x il) 

[ uf+ < uo,</?(0, .) >-s.s^ 0, 
/o Jn 

where (p is the weak solution (see Remark 4) of 



ipt + dnifi + sAip = f in (0, T) x fl, 

e{ipt ~ d,ip) - if = on(0,T)xFo, 

ipt-d^ip^O on(0,T)xFi, 

(p{T, .) = in n. 
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Using the Ricsz representation theorem and continuity (see Remark 4) of 

/ e L\{o, T)xn)^ cp{o, .) e H\n), 

it is obvious that, for any uq G X~^, there exists a solution by transposition u G i^((0, T) x 57) such that 

ll"llL2((o,T)xn) ^ IIi^oIIa'— • 
Moreover, if mq G X = Lemma [5] shows that the solution u satisfies 

l"llL2((0,T),ffi(n)) ^ ll'"o||jfO- 

If one uses classical interpolation results (see [7]), we can deduce that for every 9 g] 1/2,1], for any 
Wo G X"^^ = there exists a solution u G L'^{{0,T), H''{n)) such that 

l"llL2((o,T),if''(fi)) ^ ll^^ollA'-'d-f)- 

Using classical trace result, we have that 

\\u\\Lmo,T)xdn) < llwollx-^d-e)- (8) 

On the other hand, estimate ^ imply in particular that 

\\ua\\x <\\u\\L2{[o,T)xri) Vuq G P(f7). (9) 
Finally, ([5]) and © yields the contradictory inclusion X^''^^^^' ^ X. D 

2 Carleman inequality in dimension 1 

In this paragraph, we will establish a Carleman-type inequality keeping track of the explicit dependence 
of all the constants with respect to T and e. As in [J, we introduce the following weight functions: 

ri{x) := 2L + X, a{t, x) — — ■, (t>{t,x) 



eH{T~tY ' - eH{T-ty 

where C > e^^ . 

The rest of this paragraph will be dedicated to the proof of the following inequality: 

Theorem 9 There exists C > and sq > such that for every e G (0, 1) and every s ^ SQ{eT + e^T^) 
the following inequality is satisfied for every ipx G X: 



(10) 



J(0,T)x(-Lfi) J(0,T)x(--L.O) J(0,T)x{0,-L} 

J{0,T)x{0} 

Here, (p stands for the solution of (S") associated to (fx- 

Remark 6 One can in fact obtain the following Carleman estimate with control term in Ti 
[ q^'e-^n^\' + s^ [ ^e-^n^,f + s^[ 

J(0,T)x(-L,0) J{0,T)x{-L,0) J (0,T) x {0,-L} 

J{0,T}x{-L} 

simply by choosing the weight function i]{x) equal to x ^ — a; + L. The proof is very similar and the 
sequel will show how to deal with this case too. 
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In order to perform a Carlcman inequality for system (S"), we first do a scaling in time. Introducing 

T :~ sT, Lp{t,x) :~ ip{t/e^x) and the weights a{t,x) :~ a{t/e,x), (j)(t^x) :~ (j){t/e,x)] we have the 
following system: 

(pt+e~'^LPx + (pxx^^ in (?, 

e'^{(pt - e^'^du<p>) - (p ^0 on (7o, qj^n 

(pt — e^^dyip = on (Ti, 

(^\t=T'^V>T in(-L,0). 

if q := {0,f) X {-L,0) , a := dil x (-L,0), ctq To x (-i,0) and tri := Tj x (-L,0). We will now 
explain how to get the following result. 

Proposition 10 There exists C > 0, such that for every e S (0, 1) and every s ^ C(T + + e^^T^ + 
^.1/32^2/3^^ fliQ Jollowing inequality is satisfied for every solution of Ul\) associated to (pT G X: 



(12) 



Observe that Proposition [TO] directly implies Theorem [^1 

All along the proof we will need several properties of the weight functions: 



Lemma 11 

• ax = ~(t>, axx = ~(i>- 

There are essentially two steps in this proof: the first one consists of doing a Carlcman estimate very 
similar to that of [1]; in the second one, we will study the boundary terms appearing in the right hand 
side due to the boundary conditions. We will perform the proof of this theorem for smooth solutions, so 
that the general proof follows from a density argument. 

We will first estimate the left hand side terms of like in the classical Carleman estimate (see [3]). 
We obtain: 



Proposition 12 There exists C > 0, such that for every e G (0, 1) and every s ^ C{T + + e ^T^ + 

£i/3j'2/3-j^ i/ie following inequality is satisfied for every solution of Ul\) associated to (pT G X: 

s' I -fe-^^'^m-'+s^ 1 0e-2^"|^.p + ,s-i / r^e-2^"(|^..p + |^*n 

J, A J, ^^3^ 



Proof 

Let us introduce ip := (^e~^". We state the equations satisfied by ^. In (0,T) x (— L,0), we have the 
identity 

Pl^ + P2^ = P3^, 

where 

Flip ^ipt + 2saxipx + £"Vx, (14) 
-P2V' = ''Pxx + s'^al'ip + sat/tjj + e^^sax-tp, (15) 

and 

P31P = sa.xx'ip- 
On the other hand, the boundary conditions are: 

e^iipt + sat-ip - £~'''{'ipx + saxip)) - ip = on .x = 0, (16) 
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eijpt + satip) + -ipx + sa^ip — on a; = — L. (17) 
We take the norm in both sides of the identity in q: 

WPli'WhM + ||P2V^|[i2(,) +2(PiV,^2V')l=(,) = Il^3^lli2(,). (18) 

Using Proposition 111! we directly obtain 

Jq 

Wc focus on the expression of the double product (Piip , P2'4') {q) ■ This product contains 12 terms 
which will be denoted by Tij{ip) for 1 ^ i ^ 3, 1 ^ j ^ 4. We study them successively. 

• An integration by parts in space and then in time shows that, since ?A(0, .) = ipiTi •) = 0: '^s have 



J (Jo 

Now, we use the boundary condition (fTC|) . Wc have 



Thanks to the fact that ■0(0, .) = ■) = 0, the same computations can be done on x = —L using 
pT]) so we obtain the following for this term: 

J a J a 

Integrating by parts in time, we find 

Ti2(^) = ^ / = -s' / 5,a,t|^AP, 



and using Proposition II 11 we get 



Jq 

Tiiii') > -s-'fs 1 4>''\4>?- 

Jq 

• Now, integrating by parts in space, we have 

T2i{ip) = s / ax{\'tpx\'^)x = / + s {ax\ipx\'^ -si axli'xl^ 



Thanks to Proposition [TT] (the choice of r) is important here), the last term is positive. Using the 
boundary condition (|16p and Proposition II 11 we finally get 



Ja-o 
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An integration by parts in space provides 



T22W = / «^(IV'P). = -3^3 / ala.M^ + / allif 

J q J q ^ (7q 



This readily yields 



• Again an integration by parts in space gives 



J q J q J (70 

atax\'4'Y 



I at 

Using Proposition [m we obtain 



T2s{i^) > -s'f / - S^f / <^3|^|2 _ ^2f / ^3|^|2 



• The last integral concerning the second term in the expression of PitJj is 

T24W = -e-'s'- f alM%. 

Jq 

After an integration by parts in space, we get 

• We consider now the third term in the expression of Piip. First, we have 
Using the boundary condition we get 

T3lW>-e/ \^t\'~es'f'f 04|^|2_^-1^2 /• ^2|^|2 

J (To 

• Now, we integrate by parts with respect to x and we have 

• Then, using another integration by parts in x we obtain 

^33^ = / > -e-'sf f 02|^|2 _ f 02|^|2 _ ,~lf, I ^2|^|2 

J q J q (70 J (7\ 
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• Finally, arguing as before, we find 

Putting together all the terms and combining the resulting inequality with (jlSp . we obtain 

ll^'i'/'lli^c,) + + + ^2(^.0 < Mi') + J2(^.) + Mit) + L{i) (19) 

where the main terms are 

J q J (7\ J q 

the right hand side terms are 

Jl(^) = (s2 + sf2+e-lfs) I 02|^|2 + (52y^^^2) j ^3|^|2 

+ sf 2 /" ^■^\^\^ + [s'f + esf^) I <^3|^|2 + (sf + / ^2|^|2^ 

and 

J2(^.)= / l^.l', J3(Vt) = (5 + e + £'s) / i>\it\\ 

and the control terms are 



3l/,|2 

'■P 

(To J <JQ 

J CTq J (To J (To 

Let us now see that we can absorb some right hand side terms with the help of the parameter s. 

• First, we see that the distributed terms in Ji{ip) can be absorbed by the first term in the definition 
of /i(V') for a choice of s > f + (1 + e-i/2)f'3/2. 

• Then, we use the second term of Ii {ip) in order to absorb the integrals in the second line of the 
definition of Ji{i). We find that this can be done as long as s > f + ^2 + (e-1/2 + i)f^/2_ 

• Next, J2{4>x) can be absorbed by l2(ipx) by just taking s > T^. 

Then, we observe that all the control terms can be bounded in the following way: 

J (To 

as long as 

S>f{l + e + f 1/3^1/3 + f3/4^-l/4) ^ ^2(^ ^ ^-3/5 ^ ^-1/2 ^ ^-1/3)^ 

Since e < 1, it suffices to take s>f{l + e'^f). 

Next, we use the expression of Pi'0 and P2'0 (see (fT^ - (fTB|) ) in order to obtain some estimates for the 
terms ipt and ipxx respectively: 

s-^ f rVtl'<5' / 4>\ix\' + s-'e-' f r'li^xi^ + f rViV'P 

J q J q 

J q J q 



for any s > (1 + e~'^)T^ and 



Jo J-L Jq Jq 

+ [ 4>''\P2^J\'' < /03|^|2^^2 /■^|^^|2+ /|P2^|2, 



for s > f + £-if2. 

Combining all this with wc obtain 



/ + |^^|2) + ^-1 / 0-l(|^^^|2 + |^^|2) + ^3 / ^3|^|2 

g <J q 'J a I 

f 4>^\il;\^ + {s + e + eh) j 

for s > f{l + e-^f). 

Finally, we come back to our variable (p. We first remark that il'x = &^^°'{v>x + s(j)(p) and so 

J q J q J q 

Then, we have that ij^t = e~^'^{(pt — sat(p)^ hence 

J q J q J q 



for s >T. Analogously, one can prove that 

^^-'e-2^"|^,,|2<S-^ / / ^3|^|2+,2 , ^|^^|2^ 



s 



Jq 



for s > f2. 

Wc combine this with (j20p and we obtain the required result 

3 [ ^3^-2.a|^|2^^2 j ~^^-2s^\g,^\2 ^ ^^-1 j 1 £-2^'^ (| 1 2 + 1 2 ) + ,3 j ^3g-2.a|^|2 
J q J q J (7i 



/ ^5^-2-|^|2 + + / 0e-2-|^,p, 

J (To J (To 

for s > f (1 + e^^r) + ei/3f2/3 g^^^j ^g^^g ^j^g^^ 



which is true for s > T + ei/^T^/^ (recall that e < 1). 

With this result, we will now finish the proof of Proposition [TUl 

Estimate of the boundary term In this paragraph we will estimate the boundary term 



{s + s) cjye-^n^tl' 

J ao 
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First, wc observe that this integral can be estimated by 
for ^ ^ {T + T'^)- After an integration by parts in time, we get 

[ 02,-2.<i|^^|2^ 1 2 /■ (^2^-2.<i)^^|^|2_^2 / ^2^-2sa~~^^ 



2 

J <7q Co 



(22) 



In order to estimate the second time derivative at x = 0, we will apply some a priori estimates for the 
adjoint system. Indeed, let us consider the following function 

Then, this function fulfills the following system: 

C,t+e-\^+C,.,=em in (0,f) X (-L,0), 

- e-^dX) - C = e^Om on (0,f ) x {0}, 

Ct-e~^d^C^Ot(pt on (O,!-) X {-L}, 

C|t=T = in(-L,0). 

Using Remark |4] for (t, x) i — > Ci^t-i ^) "^c find in particular 

e [ \Ct\'< f\9t\'m'+e[i9tm\'. 

J (J J a ^ a 



This directly implies that 



Integrating by parts in time in the last integral, we have 

From the definition of 6(t) and multiplying the previous inequality by s~^, we find that 
Using Cauchy-Schwarz inequality in the last term of the right hand side of (j22p . we obtain 

J (To cr J <tq 

Combining this with ([23|) and (|2T|) yields the desired inequality (fT2l) . □ 

3 Observability and control 

In this section, we prove Theorem [TJ 
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3.1 Dissipation and observability result 

Our first goal will be to get some dissipation result. 

Proposition 13 For every e € (0, 1), for every time ti, t2 > such that t2 — ti > L and for every weak 
solution (p of (S"'), the following estimate holds true 

Mt,)\\x ^ exp |_i^l-^-^L_^| Mt2)\\x. 

Proof We first consider a weight function p{t, x) = exp(-x„) with r G (0, 1) some constant which will 
fixed later. We will first treat the strong solutions case and, using a density argument, we will get the 
weak solutions case. 

We multiply the equation satisfied by ip by pip and we integrate on Q. We get the following identity: 



We then integrate by parts in space, which due to Vp = ^/3e„, provides 
and 

-e / pipAip = e / p|V(^p + ^ / pdn{\<p\^) - e [ / pipdnf - / pipdnf 

= £j^^p\S7ip\''^L.J^p\^\'^ + L(^J^ pW?- p\p\^^-eQ'^ ppdnV-J^ ppdnV^ 
Using now the boundary conditions for p and summing up these identities, we finally get 

i ( I pW?) > '^^^ I pW\' + (1 - r) I pW\' -2e f pip,^. 



dt 

On the other hand, it is straightforward that 



d_ 

dt 



and, consequently, using that ?' G (0, 1), we have obtained 

-(1 -r) 



I (lIv^^WII^) > '-^^W^/KTMtm- 

Gronwall's lemma combined with exp(— |i) ^ p(-) 1 successively gives 

W^MWx ^ cxp - 1,)) 117^0^(^1)111 



and 

< exp (^-i (r(l - r){t, - t,) - rL)^ Mt^Wx 

We finally choose 

to — ti — L , , 

which gives the result. □ 
We will now use this dissipation estimate with our Carleman inequality to get the desired result. 
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Proposition 14 If n = 1 and for — 1 sufficiently large and e sufficiently small, the observability 
constant Ct{£) is bounded than 

C exp 



e 

where C, k are some positive constants. 
Proof 

• We begin by estimating botli sides of the Carleman inequality obtained above. We use the same 
notations as above and we define m = C — e^^ and M = C — . We first get 



T 







1 

isf) 

for tlic left liand-sidc of the Carleman inequality ([T2 



On the other hand, using that > - — — on -^j , we have the following estimate from below 



Consequently we get that 

MmUt < s%eT)-^^ exp ^^^^ 0) 0). 

We now choose s ^ (e?")^ + {sT). The above constant G is consequently estimated by 

for c' > c and c well-chosen. 
• We now deduce the result using dissipation estimates. We have just proven 



Jo JFo 

We now use the dissipation property with ti = and t2 = t E ]^,^[. We easily get, provided 
T > AL, 

which gives the resuh with = - fl - — j (T - 4L) - c' > provided that — > 8 + 32c'. 

□ 



3.2 Proof of Theorem 1 

To show the controllability result for {Sy), we will adopt some minimization strategy inspired by the 
classical heat equation. 

Proposition 15 A necessary and sufficient condition for the solution of problem (Sy) to satisfy u{T) = 
is given by: 

\/tpT <^ X,< ip{0),Uo >x= / fv. 
where ip is the solution of problem (S') with final value ipT- 
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Proof Wc apply the definition of solution to u against ip strong solutions of (S"). which gives 



ipv ^ e (piO)uo ~ s (pTu{T) + / <^(0)iio - / ^tu{T) 
"To Jdn Jon Jn Jn 

and we get the desired equivalence by approximation of weak solutions by strong solutions. □ 

Theorem 16 The following properties are equivalent 

3Ci > 0,yipT £ X; \\ip{0)\\x < Ci ||'i5||L2(-(o.T)xro) where (p is the solution of problem {S'), 



36*2 > 0,\/uq E X,3v £ L ((0, r) X Fq) such that 1|w||l2((q t^-jxTo) ^ C2||mo1|j'(: md the solution u of 



problem (Sy)satisfies u{T) ~ 0. 
Moreover, Ci = C2. 

Proof (^) Let uq G X. We define H as the closure of X for the norm defined by 



T 

2, 



\Wt\\h= ^ W{t,x)\'dadtj , 

where ip is the corresponding solution of (S"). Using the observability assumption and backward unique- 
ness (Proposition [8]), one sees that it is indeed a norm on X. 
We define a functional J in the following way 

Ji^Pr) = I I 'p'^{t,x)da-dt — f uo{x)(p{0,x)dx ~ e j uo{x)(p{0, x)da. 



2 Jo 



J is clearly convex and our assumption imply that J is continuous on H. Moreover, thanks to our 
observability assumption, J is coercive. Indeed, one has: 

for ifT G H. 

Thus J possesses a global minimum ipT G H, which gives, writing Euler-Lagrange equations, 

V^TGiJ, / / [ un^{0)+e [ uoip{0). (24) 

Jo JFo Jn Jon 

According to Proposition 1151 we have shown the existence of an admissible control defined by v = 
<^|(o,T)xro- Moreover, choosing ipx = <Pt in (|24p . we obtain the following estimate 

Hi„„ « (14' (^»s)* {i^i^>y 

Using our hypothesis, we are done. 

(<;=) If V is an admissible control with continuous dependence on uq, Proposition 1151 gives us, for every 
ipT € X, ^ 

/ uoip{0)+e / uo<^(0) = / / (pv. 
Jn Jdn Jo JFo 

Choosing now uo = (p{0) gives us the estimate 

y{Q)\\j, ^ C2yT\\H\\no\\x 

that is 

MO)\\x ^C^WiprWH- 

□ 
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Conclusion and open problems 

The question of controllability in higher dimension is open and seems to be hard to obtain using Carleman 
estimates. Another approach to obtain the observability estimate in the n-dimensional case could be 
similar to the one of Miller (see [S]) but it seems that our problem is not adapted to that framework (one 
can check that the link between the one-dimensional problem and the general one is not as simple as it 
may seem). 

If n = 1, one can wonder what is the minimal time to get a vanishing control cost when the viscosity 
goes to zero. The intuitive result would be L, but it seems that Carleman estimates can not give such a 
result. 
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